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FOREWORD 


The  work  performed  under  Contract  DAAK40-78-C-0031  that  resulted  in 
specific  deliverables  (principally  computer  programs)  is  described  in  two 
volumes  of  a  final  report.  This  is  Volume  1,  which  describes  the  work  that 
led  to  the  distributed  clutter  model.  Volume  2,  entitled  "Extended  Target 
Model:  Real-Time  Simulation  Program,"  describes  the  work  relevant  to  that 
subject. 

The  AP120B  array  processor  programs  described  in  this  report  (but  docu¬ 
mented  elsewhere)  were  designed  and  developed  by  Dr.  I.  P.  Bottlik,  as  well 
as  the  Interfaces  to  the  host  computer  and  Digital  Signal  Generator. 

Dr.  Bottlik  was  assisted  by  Mr.  D.  L.  Brandon. 
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1.  INTRODUCTION  AND  SUMMARY 

This  report  describes  a  real-time  simulation  procedure  for  generating 
replicas  of  the  gound  clutter  signal  as  received  on  the  three  monopulse  chan¬ 
nels  of  an  airborne  radar.  The  procedure  is  based  on  first  computing  the 
power  spectra  of  the  clutter  (on  the  three  monopulse  channels)  according  to 
the  engagement  geometry,  then  randomizing  these  power  spectra,  and  finally 
executing  an  FFT  in  real  time  to  create  the  quadrature  video  components  of 
the  clutter  sequences. 

^The  procedure  permits  any  spectral  shape  to  be  simulated,  so  that  it  is 
Inherently  very  accurate.  However,  the  computation  is  time-consuming  so  that 
its  application  in  real  time  is  limited.  With  one  AP120B  array  processor  in¬ 
terfaced  to  a  Datacraft/6  computer,  only  a  single  range  gate  of  clutter  can 
be  simulated  in  real  time.^ 

>  Listings  of  a  Fortran  program  are  included.  The  program  has  been  tailored 
to  a  specific  missile  system  under  test  at  the  RFSS. 

K. 
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2.  GENERATING  CLUTTER  SEQUENCES  IN  REAL  TIME 

The  use  of  an  FFT  to  generate  finite-length  sequences  of  correlated  ran¬ 
dom  samples  is  described  in  Reference  1.  First,  a  sequence  of  independent, 
complex  random  samples  is  generated  throughout  one  frequency  repetition  in¬ 
terval.  Next,  given  an  amplitude  spectrum  of  interest  (the  square-root  of 
the  ensemble-averaged  power  spectrum)  we  multiply  the  random  samples  by  the 
corresponding  samples  of  the  amplitude  spectrum.  Finally,  the  complex  spec¬ 
tral  sequence  is  Fourier  transformed  (via  an  FFT)  to  obtain  the  quadrature 
components  of  the  correlated  time  sequence. 

The  above  time  sequences  have  the  property  that  they  are  circularly 
periodic;  i.e.,  the  sequence  repeats  with  a  period  1/Af  where  Af  is  the  sam¬ 
ple  spacing  in  the  frequency  domain.  If  the  random  process  is  stationary 
we  can  make  Af  small  enough  (at  least  in  theory)  so  that  the  repetition 
period  will  be  longer  than  the  process  time  of  interest.  However,  if  the 
process  is  nonstationary,  as  with  clutter  received  on  a  missile  radar,  the 
power  or  amplitude  spectrum  of  the  process  will  be  continually  changing  and 
any  procedure  based  on  a  single  FFT  will  not  lead  to  the  proper  time  sequence. 

One  procedure  that  does  work,  however,  is  to  take  a  sequence  of  FFT's 
where  the  spectrum  for  any  one  FFT  will  be  relatively  stationary.  The  pro¬ 
cedure  is  discussed  in  Reference  2  (Section  6) .  In  order  to  avoid  the  abrupt 
transition  that  would  result  if  one  resulting  time  sequence  directly  followed 
another,  it  is  imperative  that  the  sequences  overlap  as  shown  in  Figure  1. 
Different  weighting  functions  are  discussed  in  Reference  2,  and  the  conclu¬ 
sion  is  that  linear  weighting,  as  shown  in  Figure  1,  is  satisfactory  with  an 
overlap  that  is  about  1/8  of  the  FFT  length.  More  overlap  would  provide 
slightly  lower  residual  sidelobes  but  at  the  expense  of  a  more  time-consuming 
computation. 


overlap  — *-| 


sequence  length] 
for  one  FFT 


time 


Figure  1.  Sequence  of  Weighting  Functions 


3.  GROUND  CLUTTER  SPECTRA 


The  ground  clutter  spectra  on  the  three  monopulse  channels  are  deter¬ 
mined  by  the  engagement  geometry  and  the  antenna  pattern  shapes.  In  this 
section  we  will  derive  the  basic  relationships  for  a  pulse-Doppler  radar 
with  the  antenna  pointing  within  about  30°  of  the  velocity  vector. 

3.1  Engagement  Geometry 

• 

In  Reference  2  (Section  4)  it  was  shown  that  the  clutter  spectral  shapes 
on  the  three  monopulse  channels  are  functions  of  only  five  angles,  where  the 
geometry  is  given  in  Figure  2: 

1)  azimuth  angle  of  antenna  boresight,  a 

2)  depression  angle  of  antenna  boresight,  6 

3)  depression  angle  of  range  ring  of  interest,  y 

4)  missile  dive  angle,  5 

5)  missile  roll  angle,  ui 

Two  additional  parameters  are  required  to  complete  the  specification  of  the 
clutter  spectra: 

6)  missile  altitude,  h 

7)  missile  velocity,  V 

The  engagement  geometry  may  not  be  expressed  directly  in  terms  of  the 

five  angles  above.  Instead  of  (a, 8)  it  is  most  likely  that  the  azimuth  and 

elevation  of  the  antenna  boresight  (0Q,  4>0)  relative  to  the  missile  velocity 

vector  will  be  given.  The  transformation  of  (0  ,  <b  )  to  (a, 8)  is  derived  in 

o  o 

Appendix  A. 

The  missile  altitude  and  slant  range  are  related  by 

r  -  h/sin  8  (1) 

3.2  Assumed  Antenna  Patterns 

In  the  analysis  that  follows  we  will  assume  the  following  set  of  re¬ 
ceive  voltage  gain  patterns: 

sum  channel:  g^CG ,<}>)  -  gj.(8)  (2) 

AAz  channel:  8^(9 »$)  *  g&(9) 


(3) 


Figure  2.  Engagement  Geometry  for  Ground  Clutter 
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AE1  channel:  g^E(0,4>)  -  gj.(0)  g^(<l>)  W 

where  8  and  <p  are  the  azimuth  and  elevation  angles,  respectively,  and 


gj.CO  "  exp  (-a  j  £  |  ; 


g^CO  -  kC  exp  (-c | C 1  ) 


(5) 

(6) 


where  £  designates  either  angle.  We  also  asstime  the  specific  parameters 


k  »  .084 

^  a  -  .00389 

b  -  2.30  A 
c  -  .00309 

d  -  2.10 

#  where  C  is  in  degrees  in  (5)  and  (6) .  The  half-power  width  of  the  sum  pat¬ 

tern  is  about  14°.  It  is  assumed  that  the  radar  transmits  on  the  sum  channel. 
The  patterns  for  the  sum  and  one  difference  channel  are  plotted  in  Figure  3. 

3.3  Spectral  Shape 

In  Reference  2  (Section  4)  it  was  shown  that  the  spectral  shapes  of  clut¬ 
ter  are  primarily  a  function  of  two  parameters,  ci  and  uj.  The  remaining  param¬ 
eters  in  the  list  in  Section  3.1  act  as  scaling  parameters,  either  in  frequency 
a  or  power . 


Note: 

Throughout  the  remainder  of  this  report 
clutter  spectra  will  be  shown  for  the 
three  monopulse  channels.  Unless 
otherwise  indicated,  all  plots  will 
be  in  the  format  of  Figure  4. 


/ 


Figure  3.  Assumed  Sum  and  Difference  Channel  Antenna  Patterns 


Figure  4.  Format  for  Plots  of  Clutter  Spectra 


The  azimuth  angle,  a,  strongly  affects  the  spectra  on  all  three  mono¬ 
pulse  channels  as  is  shown  in  Figure  5  through  8  for  a  »  0,  10°,  20°,  and 
30®,  respectively  (In  all  cases  g  =  5  »  uj  ■  0  and  y  *  10®,  and  frequency 
is  relative  to  2V/X).  Note  the  different  scales  in  frequency.  When  the 
azimuth  angle  is  small  (less  than  about  50%  of  the  antenna  half-power  beam- 
width)  ,  all  spectra  are  truncated  sharply  at  the  relative  Doppler  frequency 
given  by  cos  (S-y) ,  which  in  the  case  of  the  figures  is  0.985.  It  is  this 
strong  variation  with  respect  to  the  azimuth  angle  that  greatly  limits  the 
use  of  a  single  spectral  shape  to  describe  the  clutter  spectra. 

The  missile  roll  angle,  a),  affects  primarily  the  monopulse  difference 
channels  as  is  shown  in  Figures  9  through  17  for  -45®  <  w  <  45®  in  steps  of 
11.25®  (in  all  cases  8  =  <5  and  a  ■  y  *  10®).  Due  to  the  symmetry  in  the  an¬ 
tenna  patterns  of  Section  3.2,  the  spectra  for  u  and  to  +  180°  are  identical. 
The  spectra  for  (o  and  u>  +  90®  are  also  identical  except  the  monopulse  dif¬ 
ference  channels  are  interchanged  (e.g.,  see  Figures  9  and  17).  Note  that 
the  sum  channel  spectrum  is  not  a  strong  function  of  the  missile  roll  angle. 

3.4  Frequency  Scaling 

The  maximum  Doppler  frequency  within  the  range  gate  of  interest  is  given 
by 

'mx  *  ¥■  cos  W-Y>  (7) 

In  Reference  2  (Section  4.6)  it  was  suggested  that  (7)  represents  a  valid  fre¬ 
quency  scaling  of  the  spectra  for  most  cases  of  interest.  In  other  words, 
spectra  computed  at  one  combination  of  6  and  y  could  be  simply  scaled  in  fre¬ 
quency  for  other  values  by  the  use  of  (7).  The  approximation  is  good  as  long 
as  5  and  y  are  not  large. 

In  order  to  demonstrate  how  well  the  spectra  can  be  scaled  we  show  in 
Figures  18  through  23  cases  for  which  ( <5— y |  *  10®  (a  ■  10®,  8*6.  and 
<o  -  0).  The  first  three  cases  where  6  <  20®  and  y  <  30®  are  practically 
identical.  In  Figure  21  there  is  some  deviation  in  the  tail  of  the  AA  spec¬ 
trum,  but  otherwise  the  spectra  are  still  within  a  few  dB  of  the  other  cases. 
It  is  only  when  we  go  to  the  steeper  geometries  in  Figures  22  and  23  that 


A 
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Clutter  Spectra  for  a  «  20°  Figure  8.  Clutter  Spectra 


Clutter  Spectra  for  w  ■  22.5°  Figure  16.  Clutter  Spectra 


Figure  17.  Clutter  Spectra  for 


Figure  20.  Clutter  Spectra  for 
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Che  differences  become  significant.  These  results  indicate  that  the  simple 
frequency  scaling  law  of  cos  (<5-y)  applies  with  reasonable  accuracy  as  long 
as  6  and  y  do  not  exceed  about  35  or  40°. 

3.5  Clutter  Power  for  Sum  Channel 


With  the  use  of  Figure  24,  the  power  received  from  the  clutter  cell  of 


area  r  Ar  Act  is  given  by 
S  © 


AP. 


R 


PTG2x2goCrgArgAct) 
(4ir)3  r4 


(8) 


where  PT  is  the  peak  transmit  power,  G  is  the  one-way  power  gain  in  the  direc¬ 
tion  of  the  clutter  cell,  A  is  the  wavelength,  r  is  the  range,  and  0Q  is  the 
clutter  backscatter  coefficient  (radar  cross  section  per  unit  area  on  the 

ground).  The  ground  range,  r  ,  is  related  to  the  slant  range,  r,  by 

R 


g 


r  cos  y 


(9) 


The  intersection  of  the  range  resolution  cell  on  the  ground  is  Ar  .  It  is 

O 

straightforward  to  show  that 


r  Ar  *  r  Ar 
g  g 


(10) 


where  we  define  Ar  to  be  the  range  resolution  (range  gate  width). 
(8)  becomes 


Thus 


AP 


R 


P_G2\2a  ArAa 
T _ o _ 

3  3 

<4ir)J  r 


(11) 


The  total  received  power  is  obtained  by  integrating  (11)  over  all  azimuth 
angles  within  the  range  gate  intersection  on  the  ground.  If  Oq  is  constant 
(as  it  will  be  for  homogeneous  clutter) ,  the  total  power  is 


PR- 


P-A  O  Ar 
X  o 

3  3 

(4*)  r 


/ 


G  (a)  da 


(12) 


Figure  24.  Clutter  Geometry 


where  G(a)  is  defined  as  the  one-way  power  gain  in  the  direction  of  a  point 
at  an  azimuth  a  within  the  range  gate  intersection  on  the  ground. 

3.6  Effective  Radiated  Power 

In  the  RFSS  the  clutter  signal  will  be  radiated  from  the  RFSS  array  at  a 
distance  D  from  the  radar  receiver.  Let  us  designate  the  effective  radiated 
power  as  P Thus  the  power  density  at  the  receive  antenna  will  be 

Pd  -  Pe/4irD2  (13) 

The  received  power  will  be 

PR  -  PdA  -  PeA/47TD2  (14) 

where  A  is  the  antenna  aperture  area.  Since  the  peak  gain  is  given  by 
2 

G  »  4irA/A  ,  we  can  write 
t  0 


A 
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Pe-PR 


(4tt)2D2 

G  X2 
o 


(15) 


so  chat  wich  (12)  ve  obtain 


P_(J  ArD  ,  . 

p  “  •  S  /  G  <a)  dct 

4itG  rJ 
o 


/ 


(16) 


It  is  more  convenient  to  work  with  normalized  antenna  patterns,  so  we  define 


Gn(a)  -  G(a)/Go 


(17) 


where  Gq  is  the  peak  gain  on  the  sum  channel.  Thus  (16)  becomes 


(P  G D‘)  a  Ar  f  7 

p  .  — T— °-  ■  -  /  G2  (a)  da 

*  4irr3  J  n 


(18) 


which  is  our  final  result.  In  some  cases  a  is  defined  as  a  function  of  the 

o 

grazing  angle,  y,  such  as  a  sin  y  dependence. 

For  several  types  of  antenna  patterns,  it  has  been  empirically  determined 
that  for  the  sum  channel 


/ 


2  2 

G  (a)  da  *  const  *  G  (y-8) 
n  n 


(19) 


which  is  accurate  to  within  about  0.2  dB  for  all  engagement  geometries.  For 
the  antenna  pattern  defined  in  Section  3.2  the  constant  has  a  value  of  .20  so 

that 


(P  G  O  a  Ar  , 

P.  -  0.16  ■  -  3-  - -  G‘(y-B) 


(20) 


The  integral  in  (18)  is  the  quantity  "POWER"  computed  by  a  subroutine 
/(SPCTRA)  written  by  the  author. 
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Note  that  the  antenna  patterns  in  (2)  through  (6)  are  defined  in  the  normalized 
form. 

3.7  Clutter  Spectral  Parameters 

Let  us  refer  to  Figure  2.  As  viewed  from  the  radar,  the  range  ring  in¬ 
tersection  with  the  ground  plane  will  be  approximately  a  straight  line,  at 
least  for  y  not  too  large  and  for  points  on  the  range  ring  that  are  not  too 
far  away  from  the  vertical  plane  ROA.  Thus  we  have  the  approximate  geometry 
given  by  Figure  25,  where  points  A  and  G  represent  vectors  RA  and  RG  in  Figure  2, 
respectively.  The  apparent  angle  between  the  two  points  is  A  *  y  -0.  The 
skewed  coordinates  in  Figure  25  represent  the  principal  planes  of  the  antenna, 

0  being  the  azimuth  axis  and  <f>  the  elevation  axis.  The  angle  £  is  the  angle 
between  the  4>-axis  and  the  vertical  plane,  which  is  derived  in  Appendix  B  as 
a  function  of  the  engagement  geometry.  Point  Q  is  some  arbitrary  point  along 
the  ground  "line"  at  an  angular  distance  n  from  Point  G.  In  the  following 
discussion  we  will  assume  that  A,  the  angular  distance  between  Points  A  and 
G,  is  smaller  than  about  20°  so  that  the  actual  spherical  geometry  reduces  to 
the  plane  geometry  in  Figure  25. 

The  two-way  antenna  gain  pattern  for  any  of  the  three  monopulse  channels 
is  assumed  to  be  the  product  of  two  functions  (see  Section  3.2)  as 


G(9 »♦)  -  G1(6)  G2($) 


(21) 


where  G^(4>)  and  G^ (4>)  are  given  in  the  following  table: 


Channel 

sum 


A-A z 
A-El 


G^S)  G2(d>) 

gj(6)  gj(*) 

gj(9)g^(e)  gJw) 

gj(0)  g^OiOgjOJO 


with  gj.(  )  and  g^(  )  being  defined  in  (5)  and  (6),  respectively.  From 
Figure  25  we  can  write  for  an  arbitrary  point  Q  on  the  ground  within  the 
'range  ring  of  interest 


J 


A«y-B 


Ground  Line 


Figure  25.  Engagement  Geometry  as  Viewed  from  Radar 


9  -  -Asia  C  +  H  cos  £  (22) 

♦  ■  -Acos  5  -  ri  sin  £  (23) 

where  n  is  the  angular  distance  of  Point  Q  from  Point  G.  Thus  we  can  write 
the  two-way  antenna  gain  along  the  ground  line  as 

G(n)  *  G^(-Asin  £  +  n  cos  £)  G2(-Acos  £  -  n  sin  £) 

■  G^(Asin  C  -  n  cos  £)  G2(Acos  £  +  n  sin  £)  (24) 

where  symmetry  was  assumed  to  obtain  the  second  result  in  (24) . 

For  the  monopulse  sum  channel  the  two-way  pattern  will  be  nearly  circu 

larly  symmetric  about  Point  A  in  Figure  25,  assuming  the  pattern  in  (2). 

Therefore,  the  peak  gain  along  the  ground  line  occurs  at  Point  G,  and  the 

angle  from  the  boresight  axis  is  A  ■  y-B.  This  is  the  reason  why  the  effec 

2 

tlve  radiated  power  for  the  sum  channel  scales  as  G^Cy-6)  in  (20). 

For  the  monopulse  difference  channels  the  reasoning  is  similar,  except 
that  there  is  no  longer  any  dominant  circular  symmetry  about  the  boresight 
axis.  Instead,  the  symmetry  is  about  the  principal  antenna  axes.  For 
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example, Che  azimuth  difference  pattern  has  a  null  along  the  4>-axis,  so  that 
Point  N  on  the  ground  line  corresponds  to  a  null.  The  angular  separation  of 
the  null  from  Point  G  Is 


nN  -  Atan  5  (25) 

#  By  referring  to  Figure  2  we  see  that  the  azimuth  angle  of  the  null  is  a  +  ri^. 

The  normalized  Doppler  frequency  of  the  null  (ip  is  the  corresponding  Doppler 
angle)  is 


cos  *  cos  y  cos  6  cos  (a  +  n„)  +  sin  y  sin  6  (26) 

N  N 

If  we  use  the  approximation  derived  in  Reference  2  (Section  4.6),  we  can 
write  for  most  cases  of  interest 


cos  ip„  *  cos (6  -y)  cos  (a  +  tt)  (27) 

N  N 

With  the  use  of  (7)  we  can  write  the  Doppler  frequency  of  the  null  as 


N 


£„x cos  <“ +  V 


For  the  azimuth  difference  channel  there  will  be  two  points  along  the 
ground  line  in  Figure  25  for  which  the  two-way  gain  will  be  a  local  maximum. 
One  of  these  points  will  lie  to  the  left  of  Point  N  and  one  will  lie  to  the 
right.  (When  £  =  ±  90°  there  will  be  only  one  local  maximum  of  significance 
which  occurs  near  n  *  0,  and  the  null  will  occur  so  far  away  from  Point  G 
that  the  small-angle  approximation  will  no  longer  be  valid.  In  effect,  the 
clutter  spectrum  will  be  unimodal  for  the  azimuth  difference  channel  when 
£  is  in  the  neighborhood  of  ±  90°,  and  for  the  elevation  difference  channel 
when  £  =  0  and  180°  as  can  be  seen  in  Figure  13.)  There  Is  no  convenient 
way  to  find  these  peaks  other  than  to  perform  a  computer  search  for  them. 

Since  the  geometry  is  specified  by  only  two  parameters,  A  and  £,  we  have  tabu¬ 
lated  the  peak  gain  of  the  azimuth  difference  channel  that  occurs  to  the  left 
of  Point  N  in  Figure  25  as  a  function  of  these  two  angles.  The  peak  gain  is 
shown  in  Table  1  (referenced  to  the  peak  gain  of  the  sum  channel)  and  the 
location  of  the  peak,  Is  shown  in  Table  2.  Because  of  the  symmetry  we 
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Table  1.  Two-way  Gain  of  Left-Side  Peak  (dB) 


ANGLE  BETWEEN  PLANES, £  (DEG). 


A (DEG) 

C 

o 

pj 

o 

30  40  50  60  70 

30 

90 

-IS 

-61.  2 

-17 

-54.  8 

-62. 

6 

-16 

-48.  9 

-56. 

1  -65. 

6 

-15 

-43.  4 

-50. 

0  -58. 

7 

-14 

-38.  4 

-44. 

4  -52. 

3 

-63.  4 

-13 

-33.  8 

-39. 

3  -46. 

4 

-56.  4 

-12 

-29.  6 

-34. 

6  -41. 

0 

-49.  9  -63.  3 

-11 

-25.  9 

-30. 

4  -36. 

1 

-43.  9  -55.  6 

-10 

-22.  6 

-26. 

6  -31. 

7 

-38.  5  -48.  6  -65.  6 

-9 

-19.  7 

-23. 

2  -27. 

■7 

-33.  7  -42.  3  -56.  6 

-8 

-17.  2 

-20. 

2  -24. 

1 

-29.  3  -36.  7  -43.  5 

-7 

-15.  1 

-17. 

6  -21. 

0 

-25.  4  -31.  6  -41.  4  -59. 

8 

-6 

-13.  3 

-15. 

4  -IS. 

2 

-22.  0  -27.  2  -35.  2  -49. 

5 

-5 

-11.  9 

-13. 

6  -15. 

9 

-19.  0  -23.  4  -29.  7  -40. 

8  -67. 

2 

-4 

-10.  8 

-12. 

1  -13. 

9 

-16.  5  -20.  0  -25.  1  -33. 

5  -52. 

1 

-3 

-10.  0 

-10. 

9  -12. 

4 

-14.  4  -17.  2  -21.  2  -27. 

5  -40. 

1 

-2 

-9.  5 

-10. 

1  -11. 

1 

-12.  7  -14.  9  -18.  0  -22. 

6  -31. 

0 

-57. 

5 

-1 

-9.  2 

-9. 

5  -10. 

2 

-11.  3  -13.  0  -15.  3  -IS. 

7  -24. 

2 

-37. 

0 

0 

-9.  1 

-9. 

2  -9. 

6 

-10.  3  -11.5  -13.  2  -15. 

7  -19. 

2 

-25. 

2 

1 

-9.  2 

-9. 

1  -9. 

2 

-9.  6  -10.  4  -11.6  -13. 

3  -15. 

6 

-IS. 

6 

-21. 

6 

2 

-9.  5 

-9. 

2  -9. 

1 

-9.  3  -9.  7  -10.  4  -11. 

6  -13. 

0 

-14. 

7 

-15. 

8 

3 

-10.  0 

-9. 

4  -9. 

2 

-9.  1  -9.  3  -9.  7  -10. 

3  -11. 

2 

-12. 

1 

-12. 

-7 

/ 

4 

-10.  8 

-9. 

9  -9. 

4 

-9.  2  -9.  1  -9.  2  -9. 

6  -10. 

0 

-10. 

5 

-10. 

8 

5 

-11.  9 

-10. 

7  -9. 

9 

-9.  4  -9.  2  -9.  1  -9. 

2  -9. 

4 

-9. 

6 

-9. 

7 

6 

-13.  3 

-11. 

7  -10. 

7 

-9.  9  -9.  5  -9.  2  -9. 

1  -9. 

1 

-9. 

2 

-9. 

2 

7 

-15.  1 

-13. 

1  -11. 

7 

-10.  8  -10.  1  -9.  6  -9. 

4  -9. 

3 

-9. 

2 

-9. 

2 

8 

-17.  2 

-14. 

9  -13. 

2 

-11.  9  -11.0  -10.  4  -10. 

0  -9. 

7 

-9. 

6 

-9. 

6 

9 

-19.  7 

-17. 

0  -14. 

9 

-13.  4  -12.  2  -11.4  -10. 

9  -10. 

6 

-10. 

5 

-10. 

4 

10 

-22.  6 

-19. 

5  -17. 

0 

-15.  2  -13.  8  -12.  9  -12. 

2  -11. 

8 

-11. 

7 

-11. 

6 

11 

-25.  9 

-22. 

3  -19. 

5 

-17.  3  -15.  7  -14.  6  -13. 

8  -13. 

4 

-13. 

2 

-13. 

2 

12 

-29.  6 

-25. 

6  -22. 

4 

-19.  9  -18.  0  -16.  7  -15. 

8  -15. 

3 

-15. 

1 

-15. 

1 

13 

-33.  8 

-29. 

3  -25. 

7 

-22.  8  -20.  7  -19.  2  -IS. 

1  -17. 

6 

-17. 

4 

-17. 

4 

14 

-33.  4 

-33. 

4  -29. 

3 

-26.  1  -23.  7  -22.  0  -20. 

8  -20. 

2 

-20. 

0 

-20. 

0 

15 

-43.  4 

-37. 

9  -33. 

4 

-29.  8  -27.  1  -25.  1  -23. 

8  -23. 

1 

-23. 

0 

-23. 

0 

16 

-48.  9 

-42. 

9  -37. 

9 

-33.  9  -30.  8  -28.  6  -27. 

2  -26. 

4 

-26. 

2 

-26. 

4 

17 

-54.  B 

-46. 

3  -42. 

9 

-38.  4  -35.  0  -32.  5  -30. 

9  -30. 

1 

-29. 

9 

-30. 

0 

18 

-61.  2 

-54. 

2  -48. 

2 

-43.  3  -39.  5  -36.  7  -35. 

0  -34. 

1 

-33. 

9 

-34. 

1 

19 

-60. 

6  -54. 

1 

-48.  6  -44.  4  -41.  3  -39. 

4  -38. 

4 

-38. 

2 

-38. 

4 

20 

-67. 

4  -60. 

3 

-54.  3  -49.  7  -46.  3  -44. 

1  -43. 

1 

-42. 

9 

-43. 

2 

,21 

-67. 

0 

-60.  5  -55.  3  -51.  6  -49. 

3  -48. 

1 

-47. 

9 

-48. 

3 

22 

-67.  0  -61.  4  -57.  4  -54. 

S  -53. 

5 

-53. 

3 

-53. 

7 

23 

-67.  9  -63.  4  -60. 

6  -59. 

2 

-59. 

1 

-59. 

5 

24 

-65, 

3 

-65. 

2 

-65. 

/ 

22 


Table  2.  Location  of  Left-Side  Peak,  (deg) 


t 
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ANGLE  BETWEEN  PLANES, £  <  DEG ) 
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-6 
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-5 
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-7.  1  -7.  3  -7.  8  -8.  4 
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-3 
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-2 
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—6.  6 
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0 
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1 
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2 
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-6.  2  -6.  0  -5.  8  -5.  5 
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3 

-6.  5 

-6.  2 
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-4.  6 

-3.  7 
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4 
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-4.  1 
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5 

-6.  5 

-6.  1 

-5.  6  -5.  1  -4.  6  -4.  2 

-3.  7 

-2.  9 

-1.  7 

0.  0 

6 

-6.  5 

-6.  0 

-5.  5  -4.  9  -4.  3  -3.  7 

-3.  0 

-2.  4 

-1.  4 

0.  0 

7 

-6.  5 

-6.  0 
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-1.  8 

-0.  9 

0.  0 

8 

-6.  5 
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-2.  4 

-1.  5 
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-6.  5 

-5.  9 
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12 

-6.  5 

-5.  9 

-5.  3  -4.  5  -3.  6  -2.  8 

-1.  9 

-1.  0 

-0.  2 
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13 

-6.  5 

-5.  9 

-5.  3  -4.  5  -3.  7  -2.  8 

-1.  8 

-0.  9 

-0.  1 

0.  0 

14 

-6.  5 

-5.  9 

-5.  3  -4.  6  -3.  7  -2.  8 

-1.  8 

-0.  8 

-0.  1 

0.  0 

15 

-6.  5 

-6.  0 

-5.  4  -4.  6  -3.  7  -2.  8 

-1.  8 

-0.  8 

0.  0 

0.  0 

16 

-6.  5 

-6.  0 

-5.  4  -4.  7  -3.  8  -2.  8 

-1.  8 

-0.  8 

0.  0 

0.  0 

17 

-6.  5 

-6.  0 

-5.  5  -4.  8  -3.  9  -2.  9 

-1.  8 

-0.  8 

0.  1 

0.  0 

18 
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0.  1 
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19 

-6.  1 
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20 
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-5.  9  -5.  2  -4.  3  -3.  1 

-1.  8 

-0.  7 

0.  3 

0.  0 

22 

-5.  4  -4.  4  -3.  2 

-1.  9 

-0.  7 

0.  3 

0.  0 

23 

-4.  5  -3.  3 

-1.  9 

-0.  7 

0.  3 

0.  0 

24 

-0.  7 

0.  3 

0.  0 

can  write  the  location  of  the  peak  on  the  right  side  of  the  ground  line 


as 


nRCA,C)  -  -nL(-A,£)  (29) 

For  the  peaks  we  also  have  the  symmetry 

t^CA.5)  -  -nL(A,  -£)  (30) 

Note  that  nothing  changes  if  £  +  180°  is  substituted  for  £.  For  the  eleva¬ 
tion  difference  channel  we  can  use  the  above  results  by  substituting  either 
£  +  90°  or  £  -  90°  for  £. 

In  converting  the  azimuth  location  of  either  peak  to  Doppler  frequency 
we  can  write  for  any  rj 


f  "  fmax  COs(Ct  +  (31) 

Note,  however,  that  the  left  peak  will  generally  be  at  a  higher  frequency 
than  the  right  peak  (as  defined  above)  so  that  fR  <  f^  (at  least  for  a  >  0) . 
When  a  <  riL  we  get  a  foldover  or  wraparound  effect  that  results  in  the  spec¬ 
trum  being  sharply  truncated  at  f  .In  this  case  only  the  lower-frequency 

max 

peak,  namely  the  one  at  f  *  fmaT  cos (a  +  nR) ,  will  be  relatively  undistorted. 

In  the  process  of  converting  a  power  distribution  in  angle  to  a  power 
distribution  in  frequency,  we  must  also  include  the  effect  of  the  Jacobian. 
Thus  the  peak  gains  in  Table  1  (in  linear  units)  must  be  divided  by  sin(a  +  n) 
where  n  is  the  corresponding  entry  in  Table  2  (provided  a  +  £  >  0. 

3.8  Two  Examples 

Let  us  refer  to  the  parameters  of  Figure  15,  for  which  a  ■  10®,  8*0, 

Y  ■  10®,  5*0,  and  to  *  22.5®.  Immediately  we  note  that  A  »  10®  and  from 

(B-3)  we  have  £  *  22.2®.  The  parameters  that  were  discussed  in  Section  3.7 

are  summarized  in  Table  3.  The  peak  gain  for  the  sum  channel  is  defined  as 
2 

Gn(Y-8),  where  Gq(  )  is  the  one-way  power  of  the  antenna  pattern,  and  the 
Jacobian  is  defined  as  l/sin(a+r|) .  The  normalized  Doppler  is  designated 
by  f  in  Table  3.  In  the  last  column  the  "corrected  gains'1  are  referenced 
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co  Che  correcced  gain  in  Che  sum  channel.  If  we  compare  che  resulcs  of 
Table  3  wich  Figure  15  we  noce  chac  Che  locaCion  of  Che  peaks  and  nulls  is 
in  excellent  agreement,  and  the  peak  values  of  the  various  lobes  agrees  to 
within  1  dB. 

As  a  second  example  we  have  arbitrarily  selected  a  *  30° ,  B  *  5° , 
y  m  10°,  6  »  15°,  and  m  *  20°.  The  results  are  summarized  in  Table  4  along 
wich  a  plot  of  the  "exact"  sprectum.  The  agreement  is  again  excellent  as 
noted  on  Che  figure. 
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4.  APPROXIMATING  SPECTRA  IN  REAL  TIME 

As  was  concluded  in  Reference  2  (Section  4.6)  one  cannot  Implement  the 
"exact"  computation  of  the  clutter  spectra  In  real  time,  at  least  In  a  prac¬ 
tical  sense.  One  is  forced  to  use  some  approximation  technique.  In  this 
section  we  discuss  one  possibility  that  is  based  on  describing  the  clutter 
spectra  by  a  few  parameters  that  are  functions  of  the  engagement  geometry. 

4.1  Description  of  Method 

In  the  format  of  the  plots  in  Figure  4  through  23,  all  spectra  are  com¬ 
posed  of  either  a  single  lobe  or  two  adjacent  lobes  where  the  lobes  are  either 
distorted  parabolas  or  triangular-shaped  functions.  For  any  of  the  lobes  we 
have  decided  to  implement  the  parametric  description  that  is  diagramed  in 
Figure  26.  First  of  all,  a  reference  level  (base)  is  established  that  is, 
say,  40  dB  below  the  sum  channel  peak.  On  the  log-scale  of  the  figure  a 
second  level  is  established  that  is  midway  between  the  lobe  peak  and  the 
reference  base.  Five  normalized  frequencies  can  then  be  determined,  f^ 
through  for  Points  1  through  5,  respectively,  in  Figure  26.  The  sixth 
parameter  in  Figure  26  is  the  peak  height.  A  seventh  parameter,  as  shown 
in  Figure  27,  describes  the  shape  of  the  lobe.  To  reconstruct  the  spectral 
lobe  in  the  format  of  Figure  26,  a  straight  line  is  drawn  between  Points  1 
and  2.  Between  Points  2  and  3  either  a  straight  line  (TYPE  ■  0)  or  a  para¬ 
bola  with  zero  slope  at  Point  3  (TYPE  >  0)  is  drawn.  Between  Points  3  and  4 
either  a  straight  line  (TYPE  <  1)  or  a  parabola  with  zero  slope  at  Point  3 
(TYPE  ■  2)  is  drawn.  A  straight  line  is  drawn  between  Points  4  and  5. 

4.2  Application  to  Specific  System 

For  the  system  based  on  the  antenna  patterns  in  Section  3.2  we  have  de¬ 
termined  the  spectral  parameters  for  all  combinations  of  a  *  0,  10°,  20°, 

30*  and  U)  *  -45°  to  45°  in  steps  of  11.25®  (in  all  cases  0  ■  6  *  0  and 
y  •  10®).  The  results  are  shown  in  Figures  28  through  59,  where  the  a.  figure 
designates  the  "exact"  case  and  the  b.  figure  designates  the  approximate  one. 
The  agreement  is  generally  within  1  dB  throughout  the  40  dB  range  for  all 
cases.  The  parameters  are  summarized  in  Tables  5  through  8.  All  power  values 


(text  continues  on  Page  49) 
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Figure  26.  Determination  of  Six  of  the  Parameters 
that  Describe  Spectral  Lobe 


TYPE  -  0 


TYPE  -  1 


TYPE  -  2 


Figure  27.  Three  Types  of  Lobe  Shapes 


Figure  31a.  Clutter  Spectra  for  a  -  0,  w  -  -11,25°  Figure  31b.  Approximation  to  Figure  31a 


Clutter  Spectra  for  a  ■  10°,  u>  ■  -33.75°  Figure  37b.  Approximation  to  Figure  37a 


Figure  43a.  Clutter  Spectra  for  a  -  10° ,  w  -  33.75°  Figure  43b.  Approximation  to  Figure  43a 


Clutter  Spectra  for  a  -  20°,  w  -  33.75°  Figure  45b.  Approximation  to  Figure  45a 


Clutter  Spectra  for  a  -  20°,  u>  ■  -11.25°  Figure  47b.  Approximation  to  Figure  47a 


Approximation  to  Figure  49a 


Approximation  to  Figure  51a 


Clutter  Spectra  for  a  -  30°,  w  -  -33.75°  Figure  53b.  Approximation  to  Figure  53a 


Approximation  to  Figure  57a 
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Table  5.  Clutter  Spectrum  Parameters  for  Sum  Channel 

(frequencies  are  normalized  to  2V/X, 

PEAK  is  in  dB) 
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Table  8. 


Clutter  Spectrum  Parameters  for  AE1  Channel 

(frequencies  are  normalized  to  2V/X, 

PEAK  is  in  dB) 
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are  referenced  Co  Che  -40  dB  level  below  che  sum-channel  peak  (i.e.,  the  sum 
channel  peak  has  a  value  of  40  dB) .  The  azlmuch-diff erence  channel  is  char¬ 
acterized  by  two  lobes,  designated  by  "right"  and  "left"  in  the  tables, 
where  the  orientation  is  with  respect  to  the  spectrum.  In  some  cases  the 
right  lobe  is  below  the  reference  level;  the  parameters  are  set  to  zero  to 
designate  this  situation. 

In  Tables  5  through  8  a  constant  value  of  A  *  y  -  8  *  10°  was  assumed. 
In  Section  3.7  it  was  shown  that  the  spectral  parameters  are  also  a  function 
of  A.  The  dependence  is  slight  for  all  parameters  except  the  peak  power 
values  for  the  various  azimuth-and  elevation-difference  channel  lobes.  Here 
the  effect  is  strong  as  we  see  in  Table  1.  To  include  the  dependence  of  A 
on  the  peak  power  values,  the  parameters  in  Tables  6  through  8  should  be 
adjusted  by  the  values  in  Table  1  in  a  procedure  that  is  consistent  with 
the  examples  in  Section  3.8. 
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5.  SIMULATION  ARCHITECTURE 

The  computer  resources  that  are  available  at  the  RFSS  to  generate  clut¬ 
ter  signals  consist  of  a  Datacraft/1  minicomputer  that  acts  as  a  host  to  the 
AP120B  array  processor  built  by  Floating  Point  Systems.  The  input  to  the 
host  computer  consists  of  the  state  of  the  engagement  geometry  (9q,  <J>q,  r,  6, 
<i),  h,  V)  and  the  output  of  the  AP120B  is  the  set  of  three  clutter  time- 
sequences.  In  deciding  which  operations  are  to  be  assigned  to  the  two  pro¬ 
cessors,  the  following  factors  were  applied: 

1)  the  array  processor  works  best  on  arrays  where  the  number  of 
operations  on  each  sample  is  few; 

2)  the  program  construction  is  costly  for  the  array  processor 
so  that  program  should  be  fairly  stable; 

3)  modifications  are  best  made  in  the  host  computer;  and 

4)  there  should  be  relatively  little  traffic  over  the  real-time 
interface  between  the  host  computer  and  the  array  processor. 

With  these  factors  in  mind,  the  assignment  of  the  various  processing  steps  to 
generate  clutter  signals  reduces  naturally  to: 

The  host  computer 

1.  derive  and/or  transform  engagement  geometry 

2.  compute  spectral  parameters 
The  array  processor 

1.  compute  power  spectral  arrays 

2.  randomize  spectra 

3.  take  Fourier  transform  to  time  domain 

4.  weight,  overlap,  and  sum  time  sequences 

5.  normalize  time  sequences  to  reference  power 

In  order  to  minimize  the  computation  load  on  the  host  computer,  one 
should  take  advantage  of  table-lookup  procedures  to  compute  the  spectral 
parameters.  The  real-time  interface  between  the  host  computer  and  the 
array  processor  consists  of,  at  most,  28  parameters  if  the  format  of 
Section  4  is  used. 

In  Appendix  C  we  describe  and  give  listings  for  a  clutter  simulation 
program  that  conforms  to  the  above  architecture. 
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APPENDIX  A 


COORDINATE  TRANSFORMATIONS  FOR  GROUND  CLUTTER 


There  are  several  coordinate  systems  involved  for  the  ground  clutter 
as  observed  by  a  missile-borne  radar.  The  transformations  from  one  system 
to  another  are  derived  here.  We  begin  with  a  transformation  derived  in 
Reference  2  (Appendix  A),  where  one  typographical  error  is  corrected. 

Arbitrary  Point  on  Ground  to  Antenna  Coordinates 

Let  us  begin  with  the  geometry  of  Figure  60  where  R  designates  the 
location  of  the  radar  at  the  origin  of  a  rectangular  coordinate  system  ref¬ 
erenced  to  the  ground,  RV  designates  the  velocity  vector  which  is  contained 
in  the  y-z  plane  by  definition,  and  P  is  the  point  of  interest.  We  can  de¬ 
fine  the  dive  angle  of  the  missile  as  6,  the  azimuth  and  depression  angles 
of  point  P  as  C  and  y,  respectively,  and  ip  as  the  Doppler  angle.  Designating 
the  direction  cosines  a  ,  a  ,  a  on  the  x,  y,  z  axes,  we  can  write 

ax  •  cosy  sin? 

a^  »  cosy  cose  (A-l) 

a  ■  -siny 
z 

To  define  missile  coordinates  we  can  rotate  the  x-axis  clockwise  (looking  in) 
by  6  so  that  the  new  y-axis  corresponds  to  the  velocity  vector.  Now  the 
direction  cosines  in  the  new  coordinate  system  are 


Note  that  cosd 


b  or 
7 


b  ■  a  cos6  -  a  sin6 
7  7  z 

b  "a  sin6  +  a  cosd 
z  y  z 


cosij/  *  cosy  cose  cos6  +  siny  sin6 


(A-2) 


(A-3) 


If  the  missile  rotates  (rolls)  by  an  amount  u>  counterclockwise  looking  out 
t  the  velocity  vector,  then  we  can  rotate  the  previous  y-axis  by  the  corres¬ 
ponding  angle.  Defining  the  new  z-axis  to  be  the  missile  "vertical"  axis, 
we  obtain  new  direction  cosines  as 


Figure  60.  Geometry  of  Arbitrary  Point  on  Ground 


c  ■  b  cosu)  +  b  sinu 
XX  z 

c  “  b 
7  7 

c  »  -b  sinu  +  b  cosa> 
Z  X  z 


(A-4) 


We  have  now  arrived  at  the  geometry  of  Figure  61  where  0q  and  are  the  azi¬ 
muth  and  elevation  angles,  respectively,  for  the  antenna  boresight  (in  this 
derivation  we  assume  that  the  longitudinal  axis  of  the  missile  is  coincident 
with  the  velocity  vector) .  To  define  a  coordinate  system  referenced  to  the 
antenna,  we  will  first  rotate  the  x-axis  in  Figure  61  clockwise  by  0^  to 
obtain  the  new  direction  cosines  as 

d  ■  c  cosS-  -  c  sin9rt 
x  x  0  y  0 

d  •  c  sin0ft  +  c  cos9n  (A-5) 

y  x  0  y  0 

d  «  c 

f  z  z 

Mow  rotate  the  x-axis  counterclockwise  by  to  obtain  the  new  direction 
cosines  as 
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Figure  61.  Missile  Geometry 


®y  “  dyCOS^Q  +  d2Sln<t>0 

ez  ■  “dysin<f)o  +  d*cos*o 


(A-6) 


Finally,  we  can  reference  Point  P  to  the  antenna  coordinates  as  in  Figure  62 
where  the  antenna  azimuth  (9)  and  elevation  (<j>)  angles  are  given  by 


tan6  ■  e  /e 
x  y 

sin$  “  e 


(A-7) 

(A-8) 


Antenna  Boresight  (a,  8)  *♦*  (9q , 4>0) 


In  Figure  60  we  can  define  5  ■  a  and  y  ■  8  so  that  Point  P  is  the  in¬ 
tersection  of  the  antenna  boresight  vector  with  the  ground  as  shown  in 
Figure  63.  The  transformation  then  follows  (A-l)  through  (A-4)  so  that 


a^  -  cosB  sina 

a  ■  cosB  cosa 
y 

a  “  -sinfi 
z 


(A-9) 


f  and 
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and 


b  ■  a 
x  x 

b  ■  a  cos<$  -  a  sin6 

y  y  z 

b  ■  a  sin5  +  a  cos5 
z  y  z 


b  cosio  +  b  sinw 
X  z 


-b  sinw  +  b  cosw 
X  z 


These  direction  cosines  define  the  A-vector  in  Figure  61,  so  that 


Antenna  Boresight  (0q  ,  d>Q ) 


tan  9 

c 

sin  $ 

c 

(01,6) 


c  /c 
X  y 


(A-10) 


(A-ll) 


(A-12) 

(A-13) 


This  transformation  is  the  inverse  of  the  one  just  derived, 
with  the  A-vector  in  Figure  61,  which  has  direction  cosines 


We  begin 


cos$  sin0 
o  o 

COS<|>  COS0 
O  0 

sintf 


(A-14) 


z  '  o 

To  obtain  the  direction  cosines  of  the  A-vector  in  a  coordinate  system 
that  has  the  z-axis  in  a  plane  vertical  to  the  ground  we  must  rotate  the  old 
y-axis  by  oj  in  the  counterclockwise  direction  (looking  in)  so  that 


b  ■  c  cosui  -  c  sinw 
XX  z 

b  »  c 

y  y 

b  "  c  sinw  +  c  cosw 
Z  X  z 


(A- 15) 


The  missile  dive  angle  6  is  referenced  to  the  horizontal  so  we  can  de¬ 
fine  a  new  coordinate  system  referenced  to  the  ground,  that  in  Figure  63,  by 
rotating  the  old  x-axis  by  6  in  the  counterclockwise  direction  (looking  in) . 
The  new  direction  cosines  are 


A 
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a 

x 

a 

7 

a 

z 


b  cos<$  +  b  sin6 

y 

-b  sin<5  + 


The  A-vector  is  now  defined  by 


tan  a 
sin  6 


a  /  a 
x  y 

-a 


/ 


cr  n 
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APPENDIX  B 

ANGLE  BETWEEN  TWO  PLANES 


Plane  1:  Vertical  plane  containing  antenna  boresight  axis 

Plane  2:  Plane  containing  antenna  boresight  axis  and  antenna  <j)-axis 


Assume  that  the  missile  dive  angle  is  6,  the  missile  roll  angle  is  u>. 


and  the  antenna  boresight  axis  is  at  an  angle  from  the  velocity  vec¬ 


tor.  The  direction  cosines  of  the  antenna  boresight  axis  in  the  ground- 
referenced  coordinate  system  (the  A-vector  in  Figure  63)  are  given  by  the 
coordinate  transformations  in  (A-14)  through  (A-16) ,  which  we  denote  as 


(a^.a^.a^.  If  we  replace  by  <j>Q  +  90°  in  these  transformations,  we  de¬ 


fine  a  vector  that  is  coincident  with  the  antenna  4>-axis.  Let  us  designate 


the  direction  cosines  of  this  vector  as  (a^.a^.a^).  A  vector  that  is  perpen¬ 


dicular  to  the  A-vector  in  the  vertical  plane  has  the  direction  cosines 


xN 


ys 


a  a  /a 
y  z  zN 


(B-1) 


zN 


■  /a^  +  a^ 


The  two  vectors  A’  and  A^,  defined  by  their  direction  cosines  (a^,a’,a^) 


and  (a^^a^.a^) ,  respectively,  are  both  perpendicular  to  the  A-vector,  which 


is  the  intersection  of  the  two  planes  of  interest.  The  cosine  of  the  angle 
between  these  two  vectors  is  the  dot  product  A'  •  A^.  The  sign  of  the  angle 
is  given  by  the  sign  of  (A’  *  A^)  •  A,  which  can  be  rewritten  as  (A^  x  A)  •  A' 
The  direction  cosines  of  (A^  x  A)  are  (-cosa,  sina,  o) ,  so  that 


(B-2) 


(A'  x  a^)  •  A  «  (A^  x  a)  •  A'  -  -a^cosa  +  a^sina 

If  (B-2)  is  positive  the  angle  between  the  two  planes  is  positive,  and  vice 

versa. 


Instead  of  beginning  with  the  coordinate  (6  ,<j>  )  we  could  have  begun  with 

o  o 


(a, 8)  as  in  Figure  63. 
use  of  (A-9)  through  (A-13),  with  the  remaining  calculation  as  outlined  above 


The  procedure  then  would  be  to  obtain  C^0»^0)  by  the 
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The  use  of  the  above  results  is  straightforward,  but  somewhat  time  con- 
^  suming  for  real-time  applications.  We  are  thus  motivated  to  find  a  faster 

method  of  computing  the  desired  angle.  Unfortunately,  the  direct  analysis 
to  find  approximate  solutions  based  on  small  angles  is  very  cumbersome.  But 
with  a  computer  one  can  fine  through  the  process  of  trial  and  error  the  basic 
£  behavior  of  this  function  of  four  variables.  After  many  iterations  we  have 

found 

£  ■  w  +  a(8  +  (6-B)cosw)  +  ^((5-6) 2  -  a2) sin  2w  -  -j  a(S-8)w2sin22u 

#  (B-3) 

as  the  angle  between  the  two  planes  of  interest ,  where  all  angles  are  in 
radians.  The  solution  has  been  evaluated  for  0  <  a  <  30°,  -20°  <  8  <  20°, 

#  -20°  <  6  <  20° ,  and  -90°  <  cu  £  90“ .  The  error  in  the  approximation  is  shown 

\ 

in  Table  9.  In  general,  the  peak  error  is  less  than  1°,  except  when  J  6—0 1 
is  large. 


t 
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APPENDIX  C 

FORTRAN  PROGRAM  FOR  GENERATING  REAL-TIME 
CLUTTER  SEQUENCES 

This  Fortran  program  generates  continuous  time  sequences  of  ground 
clutter  that  would  be  received  on  an  airborne  radar  platform  for  a  time- 
varying  engagement  geometry.  Signals  are  generated  for  the  three  monopulse 
channels •  It  is  assumed  that  these  signals  will  be  radiated  into  specific 
positions  on  the  receive  antenna  pattern  so  that  each  signal  will  be  re¬ 
ceived  by  the  proper  channel  and  rejected  by  the  others. 

This  program  was  written  for  a  single  range  gate  of  a  pulse-Doppler 
radar,  although  it  is  easily  extended  to  multiple  range  gates  (but  not  in 
real  time  with  a  single  AP120B) .  The  organization  of  this  program  is 
sketched  in  Figure  64.  At  the  RFSS  it  is  assumed  that  the  real-time  state 
of  the  environment  will  be  set  up  in  C0MM0N/C1/  and  periodic  calls  will  be 
made  to  CLUTTR  in  the  Datacraft/6.  The  subroutines  that  actually  generate 
the  clutter  sequences  (FILL,  CURVE,  PARAB,  CONVRT,  TIMESQ,  GAUSS4,  FOURT) 
will  be  implemented  in  the  AP120B.  Their  Fortran  equivalents  are  included 
in  this  package. 

A  test  program  (TEST)  is  also  included  to  compare  results  on  different 
computers.  The  output  of  this  program  is  a  power  spectral  estimate  of  the 
clutter  sequences  generated  under  stationary  conditions. 

Note:  this  program  does  not  include  an  adjustment  for  A  on  the  mono¬ 
pulse  difference  channel  spectra  as  discussed  in  Sections  3.7  and  4.2.  The 
modification  is  minor  to  include  this  adjustment,  and  only  the  portion  of 
the  program  that  resides  in  the  host  computer  will  be  affected. 


I 
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PROGRAM  SUMMARY 


Main  or  Driver 


TEST  (for  this  example) 


Initialization 

RCDATA 

TIMSQI 

GADSSI  (GAUSS3.RAND) 
Real-time  Fortran 

CLUTTR(GAIN) 

PARMS(AINT) 


Utility 


PR 

ACCUM 

POWSPT 

XMIT 

FOURT 


Real-time  AP12QB 


i 

! 


FILL 

CURVE 

PARAB 

CONVRT 

TIMESQ 

GAUSS4 


/ 


non  fioo  o  o  o  o  o  o 


64 


PROGRAM  TEST < INPUT* OUTPUT,  TAPE5=INPUT,  TAPE6=0UTPUT> 

TEST  OF  SUBROUTINE  CLUTTR 

THE  TEST  CASE  CORRESPONDS  TO  ALPHA=20.  , BETA=0.  , GAMMA=10.  ,  DELTA=0.  ,  AND 
OMEGA—22.  5  (ALL  DEGREES). 

DIMENSION  SR  (2240,  3),  SI  (2240,  3),  S(200,  3),  W0RM600) 

COMMON  /Cl/  WL,  H,  V,  R,  FR,  DE,  OM,  THO,  PHO 
COMMON  /C2/  PTCDSG, SI GO,  DR 
COMMON  /C3/  NFFT, NOVLP, NF, A (64) 

COMMON  /C5/  PEFF, 

1  XRl  (320) ,  XII  (320),  XR2(320),  XI2(320),  XR3(320),  XI3(320> 


THE  FOLLOWING  PARAMETERS  MUST  BE  SET. 


DATA  WL,  H,  R,  FR,  V/.  02,  1736.  ,  10000.  ,  10000.  ,  500.  / 
DATA  DE,  OM,  THO,  PHO/O.  ,  -.  3927,  .  3244,  .  1313/ 

DATA  PTGDSQ,  DR,  SIGO/1.  ,  20.  ,  .  01/ 

DATA  NFFT, NOVLP /256. 32/ 

HERE  WE  INITIALIZE 

CALL  RCDATA 
CALL  TIMSQI 
CALL  GAUSS I 


C  HERE  WE  GENERATE  CONTIGUOUS  TIME  SEQUENCES 
C 

L»l 

DO  20  1-1,  10 
CALL  CLUTTR 

XMIT (NF,  XRl, SR (L,  1 ) ) 

XMIT (NF, XR2, SR(L, 2) ) 

XMIT(NF, XR3, SR(L,  3) ) 

XMIT  (NF,  XII,  SKL,  1)  ) 

XMIT  (NF,  XI 2,  SKL,  2)  ) 

XMIT (NF, XI3,  SI (L, 3) > 


IN  THE  SR  AND  SI  ARRAYS 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
L-L+NF 
20  CONTINUE 


C  HERE  WE  COMPUTE  THE  (AVERAGED)  POWER  SPECTRUM 
C 

L-l 

'  DO  30  1-1, 11 

CALL  POWSPT ( SR  ( L,  1  ) ,  SI  (L,  1  > ,  200,  200,  0.  ,WORK> 

CALL  POWSPT  ( SR  <L.  2),  SKL,  2),  200,  200,  0.  .WORK) 

CALL  POWSPT  ( SR  (L,  3),  SI  <L,  3),  200,  200,  0.  ,  WORK) 

CALL  ACCUM(  SR  (L,  1),S(1,  1),  200,  I) 

CALL  ACCUM(SR(L, 2), S(l, 2),  200, I) 
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ft 


ft 


ft 


SUBROUTINE  RCDATA 
C 

C  THIS  SUBROUTINE  READS  THE  CLUTTER  DATA  BASE  INTO  COMMON  /Zb/. 

C 

DIMENSION  I W  ( 4 ) ,  W  <  22 ) 

COMMON  /C6/  S< 180) , IS (4) , A1 (216) ,  IA1 (36),  A2< 180),  IA2<36> ,  EC216) . 
1  I E ( 36 ) 

DATA  IW/10HIS . ,  10HIA1 . ,  10HIA2 . ,  10HIE . / 

DATA  W/10HS1 . ,  10HS2 . ,  10HS3 . ,  10HS4 . . 

1  10HS5 . ,  10HA11 . ,  10HA12 . ,  10HA13 . . 

2  10HA14 . ,  1  OH  A 1 5 . ,  10HA1P . ,  10HA22 . . 

3  10HA23 . .  10HA24 . 10HA25 . 10HA2P . . 

4  10HE1 . ,  10HE2 . .  10HE3 . ,  10HE4 . 

5  10HE5 . ,  10HEP . / 


READ  102,  S 
READ  101,  IS 
ft  READ  102,  A 1 

READ  101,  IA1 
READ  102,  A2 
READ  101,  IA2 
READ  102,  E 
READ  101, IE 

PRINT  200 


Kl-1 

K2-36 

DO  10  L-l, 5 

•  PRINT  202, W<L)»  (S<K>, K=K1,  K2) 

Kl-Kl+36 
K2-K2+36 
10  CONTINUE 

p  PRINT  201,  IW(1),  IS 

Kl-1 

K2-36 

DO  20  L»6»  1 1 

PRINT  202, W<L).  < A1 (K) . K-Kl, K2> 
»  Kl-Kl+36 

K2-K2+36 
20  CONTINUE 

,  PRINT  201,  IW(2),  IA1 

ft  Kl-1 

K2-36 

00  30  L -12,  16 

PRINT  202, W<L),  (A2(K), K-Kl,  K2> 
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Kl»Kl+36 
K2-K2+36 
30  CONTINUE 

PRINT  201  * IW<3>» IA2 

Kl-1 

K2-36 

00  40  L=17>  22 

PRINT  202*  W(U)<  <E(K),  K»K1,  K2> 
Kl-Kl+36 
•  K2*K2+36 

40  CONTINUE 

PRINT  201* IW(4) *  IE 


100 

101 

102 

200 

201 

202 


RETURN 

FORMAT  (20F4.  0) 

FORMAT (3612) 

FORMAT (2 <9F4.  0,  4X) ) 

FORMAT <  1H1.  50HGR0UND  CLUTTER  DATA  BASE  READ  BY  SUBROUTINE  RCDATA/ ) 
FORMAT  < /1XA10* 3X3613) 

FORMAT ( / 1 XA1 0*  3X9F7.  3/U4X9F7.  3) ) 

END 


I 
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SUBROUTINE  TIMSQI 

THIS  SUBROUTINE  PERFORMS  INITIALIZATION  FOR  SUBROUTINE  TIMESQ. 

THE  INITIALIZATION  PARAMETERS  ARE . 

NFFT  -  SIZE  OF  FFT  < =2** INTEGER > 

NOVLP  «  NUMBER  OF  SAMPLES  THAT  THE  FFTS  OVERLAP  (NFFT/ 8  IS  OK) 

THE  ARRAYS  MUST  BE  DIMENSIONED  AS  LARGE  AS . 

XR1,  XII. ETC  NFFT+NQVLP-i 

A  NOVLP -I 

COMMON  /C3/  NFFT,  NOVLP,  NF, A (64) 

COMMON  /C5/  PEFF, 

1  XR1 (320),  XII (320),  XR2(320),  XI2(320), XR3(320),  XI3(320) 

•  NF-NFFT-NOVLP 
N 1 -NOVLP -1 
N2-NFFT+N1 

DO  10  K— 1 , N1 
A ( K ) -K/FLOAT ( NOVLP ) 

10  CONTINUE 

•  CALL  XMIT(-N2, 0.  .  XRl) 

CALL  XMIT (— N2,  0.  .  XR2) 

CALL  XMIT < -N2,  0.  ,  XR3) 

CALL  XMIT ( -N2,  0.  ,  XII) 

CALL  XMIT (-N2,  O.  ,  XI2) 

CALL  XMIT (-N2,  0.  .  X 13) 

•  RETURN 
END 


ooooouuoo 
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SUBROUTINE  GAUSS I 

THIS  SUBROUTINE  FILLS  ARRAYS  TR  AND  TI  WITH  INDEPENDENT  GAUSSIAN 
RANDOM  NUMBERS  FOR  USE  WITH  SUBROUTINE  GAUSS4. 

ARRAYS  TR  AND  TI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTR+NMAX  AND  NTI+NMAX 
RESPECTIVELY. 

IF  N  BELOW  IS  NOT  EVEN,  ADD  1  TO  THE  DIMENSION  OF  TI. 

COMMON ' /CS/  IRSET, NTR, NT I,  NMAX, KR, KI, TR  < 1253) , TI ( 1353) 

•  DATA  NTR,  NTI/997,  1097/,  NMAX/256/,  KR,KI/1,1/ 

N»NTR +NT I +2*NMAX 

DO  20  K»1,N,  2 

CALL  GAUSS3(TR<K), TR(K+1 ),  1.  ) 

20  CONTINUE 
RETURN 

•  END 


uuouuu 
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SUBROUTINE  GAUSS3<  X,  Y,  P ) 

THIS  SUBROUTINE  GENERATES  RANDOM  PAIRS  FOR  THESE  DISTRIBUTIONS . 

CALL  GAUSS3( X,  Y»  P ) _ GAUSSIAN  COMPONENTS  OF  AVERAGE  POWER  P 

CALL  RANG3<X,Y,  P> . RANDOM  PHASOR  COMPONENTS  OF  POWER  P 

DATA  I/O.  / 

C  )X<FNAR»)X(DNAR 

IF<P)  30,20.10 
10  E«SQRT(-P*ALOG<RAND<Z) > ) 

CO  TO  15 
ENTRY  RANG3 
IF(P )  30,20,12 
12  E*SQRT(P) 

15  A*RAND(Z) 

A=A+A- 1 . 

B*RAND<Z) 

B-B+B-l. 

A2®A*A 

B2«B*B 

C-A2+B2 

IF(C.  GT.  1.  )  GO  TO  15 
X»E*<A2-B2)/C 
Y«E*2.  *A*B/C 
RETURN 
20  X*0. 

Y-0. 

RETURN 
30  STOP  44 
END 


f 
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FUNCTION  RAND<  Z ) 

C  FOR  TEST  PURPOSES  ONLY. 

DIMENSION  A< 8) 

DATA  A/.  1.  .  2,  .  3.  .  4.  .  5/  .  b,  .  7,  .  8/ 
DATA  K/l/ 

RAND*A(K) 

K»K+1 

IF<K.  OT.  8)  K=1 

RETURN 

END 


/ 


SUBROUTINE  CLUTTR 


C 

C  IN  THIS  SUBROUTINE  WE  GENERATE  THE  CLUTTER  MODULATION  SIGNALS  FOR  ONE 
C  RANGE  GATE  ON  EACH  OF  THREE  MONOPULSE  CHANNELS. 

C 

C  THE  BACKSCATTER  COEFFICIENT  OF  THE  GROUND  IS  ASSUMED  TO  VARY  AS  THE 
C  SINE  OF  THE  GRAZING  ANGLE. 

C 

C  THE  REAL-TIME  PARAMETERS  ON  INPUT  ARE . 

C 

C  WL  *  WAVELENGTH 

C  H  »  PLATFORM  ALTITUDE 

C  V  -  PLATFORM  VELOCITY 

C  R  -  RANGE  TO  RANGE  GATE  OF  INTEREST  (UNAMBIGUOUS) 

C  FR  *  EQUIVALENT  SAMPLE  RATE  IN  RECEIVER  (PROCESSING  BAND) 

C  DE  *  DELTA,  THE  MISSILE  DIVE  ANGLE 

C  OM  *  OMEGA,  THE  MISSILE  ROLL  ANGLE  (CCW  IS  POSITIVE) 

C  THO  *  THETA-ZERO,  THE  AZIMUTH  ANGLE  OF  THE  ANTENNA  BORESIGHT 

C  PHO  -  PHI-ZERO,  THE  ELEVATION  ANGLE  OF  THE  ANTENNA  BORESIGHT 

C 

C  THE  NONREAL-TIME  PARAMETERS  ON  INPUT  ARE . 

C 

C  PTGDSQ  *  PRODUCT  OF  TRANSMIT  POWER,  GAIN,  SQUARE  OF  CHAMBER  LENGTH 
C  DR  -  RANGE  RESOLUTION 

C  SICO  -  GROUND  BACKSCATTER  COEFFICIENT  AT  30-DEG  GRAZING  ANGLE 

C 

C  THE  REAL-TIME  OUTPUT  FOR  THAT  PORTION  COMPUTED  IN  THE  DATACRAFT  IS 
C  THROUGH  /C4/,  AND  IT  IS  COMPUTED  IN  SUBROUTINE  PARMS. 

C 

C  SOME  OF  THE  OTHER  QUANTITIES  USED  IN  THIS  SUBROUTINE  ARE . 

C 

C  AL  »  ALPHA,  THE  AZIMUTH  ANGLE  OF  THE  BEAM  AXIS  IN  GROUND  COORD. 

C  BE  -  BETA.  THE  ELEVATION  ANGLE  OF  THE  BEAM  AXIS  IN  GROUND  COORD. 

C  CA  -  GAMMA,  THE  GRAZING  ANGLE  AT  THE  RANGE  RING  ON  THE  GROUND 

C 

C  ALL  ANGLES  ARE  IN  RADIANS  AND  ALL  DISTANCES  IN  METERS. 

C 

COMMON  /Cl/  WL,  H,  V,  R,FR,  DE,  OM,  THO.  PHO 
COMMON  /C2/  PTGDSQ,  SI GO,  DR 
COMMON  /C3/  NFFT 
DATA  FCTR/i.  / 

C 

C  COMPUTE  GEOMETRY  AND  DOPPLER  PARAMETERS 
C 

'  SGA-H/R 

CQA-SQRT< 1.  -SGA**2> 

OA“ATAN<  SGA/CGA) 

CPHO*COS  <  PHO ) 

SPHO*S I N  <  PHO ) 

CTHO*COS  <  THO ) 


uuuuou  ouu 
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STHO-S I N ( THO  > 

COM-COS < OM ) 

SOM»SIN<OM> 

CDE-COS(DE) 

SDE=SIN<DE> 

AX-CPHO*STHO 
AY»CPHO*CTHO 
AZ-SPHO 

BX»AX*COM-AZ*SOM 
BY-AY 

BZ»AX*SOM+AZ*COM 
CX*  BX 

CY-  BY*CDE+BZ*SDE 
CZ— BY*SDE+BZ*CDE 
AL-ATAN<CX/CY> 

BE*ATAN< -CZ/SQRT < 1.  -CZ**2) ) 

PEFF*.  013*PTGD5Q* ( 2.  *SIGO*SGA)*DR*GAIN(GA-BE> /R**3 
FMAX—2.  *V/WL 

FMC  =FMAX*<CGA*CDE+SGA*SDE> 

FSCL-.  9848*FMAX/FMC 
IAMB-FMC/FR 

FNP1=<  IAMB-1  )*FR*FSCL/FMAX 
DFNP-FNP1 /FLOAT <  < I AMB-1 ) *NFFT > 

FNP1-.  984S-FCTR*< .  9848-FNP1) 

DFNP»FCTR*DFNP 
ALD=57.  3*AL 
OMD-57.  3*0M 

IF(  ALD.  GT.  O.  >  GO  TO  10 
OMO— 360.  -OMD 
ALD— ALD 
OMDP-OMD 
IF(OMD.  CT.  315.  )  0MDP-0MDP-360. 

I CASE-  < OMDP+45.  ) /90. 

OMDP-OMDP - I C ASE*90. 

I ROT-O 

IF(  I  CASE.  EQ.  1.  OR.  ICASE.  EQ.  3)  IROT-1 

)/4.  21E.  4.  9F3.  1.  9F2,  1.  8F3,  4.  8F3,  413/ . RTTULC  H21/  (TAMROF 

FFEP. PNFD,  1PNF  1 

,  LCSF,  CMF. XAMF. PDMO. DMO. DLA, EB, LA. AG, TORI.  ESACI,  BMAI, 001  TNIRP 

COMPUTE  REAL-TIME  PARAMETERS 

CALL  PARMS< ALD,  OMDP. FNP 1 , DFNP ) 

GENERATE  REAL-TIME  SPECTRA 

I ROT-O 

IF<  ICASE.  EQ.  1.  OR.  ICASE.  EQ.  3)  IROT-1 
CALL  FILL(IROT) 

RETURN 

END 


10 


od 
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FUNCTION  CAIN (PHI) 

HERE  WE  COMPUTE  THE  2-WAY  POWER  GAIN  BY  TABLE  LOOKUP.  PHI  IS  THE 
ANGLE  IN  RADIANS. 

NB  »  NUMBER  OF  SAMPLES  IN  TABLE 
DB  =  SAMPLE  SPACING  ( DEG > 

THE  FIRST  SAMPLE  IS  AT  PHI=0. 

DIMENSION  BEAM < 101 ) 

DATA  NB/101/,  DB/.  5/ 

DATA  A/.  00389/,  8/2.  30/ 

DATA  NO/O/ 

IF  (NO.  GT.  0)  GO  TO  30 
N0=*1 

DO  20  K=1,NB 

BEAM ( K >  =EXP ( -4.  *A*< (K-l >*DB)**B> 

20  CONTINUE 
30  £»57.  3*ABS(PHI> 

IE-E/DB+1.  5 
GAIN=0. 

IF ( IE.  LE.  NB >  GAIN=BEAM(IE) 

RETURN 

END 


t 
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SUeROUTINE  PARMS( ALPHA,  OMEGA,  FI,  DF> 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  REAL-TIME  PARAMETERS  THAT  WILL  BE 
SENT  TO  THE  AP120B.  THE  PARAMETERS  ARE  DETERMINED  BY  INTERPOLATION 
WITHIN  THE  DATA  BASE  STORED  IN  COMMON  /C6/. 

THE  INPUT  IS . 

ALPHA  =  ABSOLUTE  VALUE  OF  AZIMUTH  ANGLE  (DEG) 

OMEGA  =  MISSILE  ROLL  ANGLE  (DEG,  -45.  .  LE.  OMEGA.  LE.  45.  ) 

FI  =  NORMALIZED  DOPPLER  OF  FIRST  SAMPLE 
DF  »  SAMPLE  SPACING  IN  DOPPLER 

THE  OUTPUT  IS  THROUGH  COMMON  /C4/ . 


IS, SPK . SUM  SPECTRUM  PARAMETERS 

IA1, A1PK.  .  .  . DEL-AZ  SPECTRUM  PARAMETERS,  LEFT  HALF 
IA2, A2PK.  .  .  .  DEL-AZ  SPECTRUM  PARAMETERS,  RIGHT  HALF 

IE, EPK . DEL-EL  SPECTRUM  PARAMETERS 

PSL . SIDELOBE  LEVEL  (SEE  SUBROUTINE  FILL) 


SKP, A1PK,  A2PK.  EPK,  AND  PSL  MUST  BE  IN  THE  SAME  UNITS  (NEPERS). 

TNIA  LANRETXE 

COMMON  /C4/  IS(6),  IA1 (6),  IA2(6).  IE(6). SPK,  A1PK,  A2PK,  EPK,  PSL 
COMMON  /C 6/  S 1 ( 36 ) ,  S2 ( 36 ) ,  S3 ( 36 ) , S4 ( 36 ) <  S5 ( 36 ) ,  I ST ( 4 ) , 

1  All (36) ,  A12( 36  > .  A13 ( 36 ) , A1 4 ( 36 ). AI5(36), A1P( 36) ,  IA1T(36>, 

2  A22  <  36 ) ,  A23  ( 36 ) ,  A24  ( 36 ) ,  A25  ( 36 ) ,  A2P  ( 36 ) »  I A2T  ( 36 ) , 

3  El ( 36 )  i  E2 ( 36 ) ,  E3 ( 36 ) , E4 ( 36 ) »  E5 ( 36 ) ,  EP ( 36 ) ,  IET(36> 

COMMON  /CINT/  K1 ,  K2,  K3,  K4,  HOM,  HAL 

DATA  NAL,  NOM/4,  9/ 

DATA  DAL,  DOM/ 10.  ,11.  25/ 

I  SAMP  ( A)  =  ( AINT  ( A )  -FI )  /DF+1 .  5 

HERE  WE  FIND  THE  PARAMETERS  FOR  INTERPOLATION 

AL -ALP HA/ DAL 

lAL-Ac 

HAL  -AL  - 1  AL 

ial»ial*i 

r  :**_  GE  NAL)  HAL-1. 
r  NA*_)  IAL-NAL-1 

DOM 


non 
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K3-K1+N0M 
K4=K2+N0M 
KAL=ALPHA/DAL+1.  5 
KAL=MINO<KAL»  NAL) 

KQM= ( OMECA+45.  )/DOM+l.  5 

•  K*(KAL— 1 ) *NOM+KOM 

HERE  WE  COMPUTE  THE  REAL-TIME  PARAMETERS 

IS(1)  *ISAMP<S1) 

IS(2)  »ISAMP(S2) 

ISO)  =ISAMP<S3) 

IS(4)  »ISAMP(S4> 

ISO)  =ISAMP(S5) 

IAi < 1 )*ISAMP (All) 

IA1(2)»ISAMP(A12) 

IAI <3)=ISAMP ( A13) 

IA1(4)»ISAMP(A14) 

IAI (5)=ISAMP  < A15) 

IA2(  1)-IA10) 

IA2<2)=ISAMP(A22) 

IA2(3)=ISAMP(A23) 

IA2(4)*ISAMP(A24) 

IA2(5)»ISAMP(A25) 

IE( 1 )  =ISAMP<E1) 

IE(2>  = I SAMP ( E2 ) 

IE(3)  * I SAMP  < E3 ) 

IE(4)  =ISAMP(E4) 

IE(5)  =ISAMP(E5) 

SPK  =»40. 

PSL  *0. 

A1PK»AINT(AIP) 

A2PK«AINT(A2P) 

EPK  =»A  I  NT  <  EP  ) 

IS(6)  *IST(KAL) 

IA1(6)»IA1T<K) 

IA2(6)*IA2T(K) 

IE<6>  -IET(K> 

C 

C  CONVERT  DB  TO  NEPERS  (AMPLITUDE) 

•  C 

SPK  1 1513*SPK 
AIPK-. 1 1513*A1PK 
A2PK*  1 1513*A2PK 
EPK  ».  1 1513*EPK 
PSL  1 1513*PSL 


RETURN 

END 


o  o  o 
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FUNCTION  A I NT (A) 

INTERPOLATION  FOR  SUBROUTINE  PARMS 

•  DIMENSION  AID- 

COMMON  /CINT/  Kl.  K2/  K3*  K4,  HQM>  HAL 
Al«(l.  -H0M>*A(K1)+H0M*A(K2> 

A2»  < 1 .  -HOM >  *A  <  K3 ) +HOM*A ( K4 ) 
AINT«(1.  -HAL ) *A1+HAL*A2 
RETURN 

•  END 


o  o  uo  o  o 
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SUBROUTINE  PR(A»  NX. NY) 

IN  THIS  SUBROUTINE  WE  PRINT  THE  2-DIMENSIONAL  ARRAY  A  WITH  THE  FORMAT 
EQUIVALENT  TO  THE  DIMENSION  A<NX.NY>. 


THE  PRINTING  IS  ON  FILE  LU,  WHICH  MUST  BE  EQUIPPED. 


20 


100 

101 

102 


DIMENSION  A < 1 ) 

COMMON  /MESGE/  MFLAG 
DATA  LU/6/ 

IF ( MFLAG.  NE.  1)  WRITE  (LU,  100) 

MFLAG=0 

DO  20  J=l, NY 

L1=(J-1)*NX+1 

L2*Ll+9 

L2MAX*J*NX 

L2=*MIN0(L2,  L2MAX ) 

DO  20  1*1, NX. 10 

IFd.LE.  1)  WRITE  (LU,  101)  V,  (A(K>,  K*L1>  L2> 
IFd.GT.  1)  WRITE  (LU.  102)  ( A(K> ,  K*L1,  L2> 

L1*L1+10 

L2*MIN0(L2+10,  L2MAX) 

CONTINUE 
RETURN 
FORMAT (1H1) 

FORMAT ( /I4, 2X10E12.  4) 

FORMAT (6X10E12.  4) 

END 


t 
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SUBROUTINE  ACCUM<X, Yi N,  IREP) 

IN  THIS  SUBROUTINE  WE  ACCUMULATE  AN  ARRAY.  AT  ALL  TIMES  ARRAY  Y 
REPRESENTS  A  RUNNING  AVERAGE. 

DIMENSION  X(l). Y(l) 

IF ( IREP.  EG.  1)  GO  TO  30 
REP*IREP 
DO  20  K=li  N 

Y(K)»( (REP-1.  >*Y<K>+X<K> >/REP 
20  CONTINUE 
RETURN 

•  30  CALL  XMIT (N* X,  Y) 

RETURN 

END 
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SUBROUTINE  POWSPT(XR,  XI,  NIN,  NOUT,  ALPHA,  W> 

C 

C  IN  THIS  SUBROUTINE  WE  COMPUTE  THE  POWER  SPECTRUM  OF  THE  COMPLEX  TIME 
C  SEQUENCE  IN  THE  ARRAY-PAIR  <XR,XI>  OF  LENGTH  NIN.  THE  POWER  SPECTRUM 
C  IS  RETURNED  IN  ARRAY  XR,  AND  IT  IS  NOW  OF  LENGTH  NOUT. 

C 

C  THE  SAMPLE  SPACING  OF  THE  POWER  SPECTRUM  IS  1/NOUT  OF  THE  REPETITION 
C  FREQUENCY. 

C 

C  A  COSINE-ON-A-PEDESTAL  WEIGHTING  IS  APPLIED  TO  THE  INPUT  SAMPLES. 

C  ALPHA  IS  THE  RATIO  OF  THE  WEIGHTING  FUNCTION  AT  THE  EDGE  TO  THE 
C  CENTER.  ALPHA5*.  03  FOR  HAMMING  AND  ALPHA5*  1.  0  FOR  UNIFORM  WEIGHTING. 

C 

C  ARRAY  W  IS  A  WORKING  ARRAY  AND  IT  MUST  BE  DIMENSIONED  AS  LARGE  AS.  .  .  . 
C 

c  NIN  IF  NOUT.  EQ.  2** INTEGER 

C  NIN+2*N0UT  IF  NOUT.  NE.  2** INTEGER 

C 

C  THE  WEIGHTS  ARE  NORMALIZED  SO  THAT  THE  SUM  IS  UNITY. 

C 

DIMENSION  XR(1), XI ( 1 ) ,  W( 1 ) 

DATA  TWOP 1/6.  2831353/ 

DATA  NO/O/ 

IF  (NO.  GT.  0)  GO  TO  25 
N0=*1 

N=*NOUT— NIN 
A»(l.  +ALPHAJ/2. 

B-(l.  -ALPHA) /2. 

CN=(NIN+l)/2. 

XN»NIN 

DO  10  K=l,NIN 

W  <  K ) =A+B*COS ( TWOP I *  <  K-CN  > / XN ) 

10  CONTINUE 
SUH**0. 

DO  15  K=l,  NIN 
SUM«SUM+W(K) 

15  CONTINUE 
WNORM»SUM 
DO  20  K=l,  NIN 
W  <  K ) *W  <  K ) / WNORM 
20  CONTINUE 
25  I F(N)  60,  35,  30 
30  CALL  XMIT  <-N»  0.  , XR (NIN+l ) ) 

CALL  XMIT(-N, 0.  , XI (NIN+l)) 

36  DO  40  K*1  * NIN 
XR(K)»XR(K)*W<K) 

XI(K)»XI<K)*W<K) 

40  CONTINUE 

CALL  FOURT <  XR,  XI,  NOUT,  1,  1, 1, W(NIN+1 > . W(NIN+N0UT+1 ) ) 

DO  50  K«l, NOUT 
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XR(K>=XR<K)**2+XI<K)**2 
50  CONTINUE 
RETURN 

60  PRINT  100,NIN,  NOUT 
STOP 

•  100  FORMAT </lX49H***ERR0R  IN  POWSPT  ,  NOUT  IS  SMALLER  THAN  NIN  ***/ 

1  /10X4HNIN*I6.  10X5HN0UT»I6) 

END 


f 


J 
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SUBROUTINE  XMIT(M» A. 3 ) 

IN  THIS  SUBROUTINE  WE  EITHER  TRANSMIT  ARRAY  A  TO  ARRAY  B  (IF  N.  GT.  O) 
OR  WE  TRANSMIT  THE  CONSTANT  A  TO  ARRAY  B  (IF  N.  LT.  0).  IN  EITHER  CASE 
T!!F  ARRAY  LENGTH  IS  IABS(N). 

THIS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEM3LY  LANGUAGE 

DIMENSION  A(l). 3(1) 

IF(N>  10,20,25 
10  NN—N 
AA*-A  ( 1 ) 

DO  15  K-l.NN 
B  <K)**AA 
15  CONTINUE 
20  RETURN 
25  DO  30  K=1,N 
B (K)*A(K) 

30  CONTINUE 
RETURN 
END 


SUBROUTINE  FILL ( IROT) 


C 

C  IN  THIS  SUBROUTINE  WE  GENERATE  THE  REAL-TIME  CLUTTER  MODULATION 
C  SIGNALS  ON  THE  THREE  MONOPULSE  CHANNELS  (FOR  ONE  RANGE  GATE  ONLY). 

C  THIS  SUBROUTINE/  AS  WELL  AS  SUBROUTINES  CURVE/  CONVRT.  PARAB,  TIMESQz 
C  GAU3S4,  AND  FFT2z  WILL  BE  IMPLEMENTED  ON  THE  AP120B. 

C 

C  THE  REAL-TIME  INPUT  IS  THROUGH  /C4/ . 


C 

C  IS,  SPK . SUM  SPECTRUM  PARAMETERS 

C  IA1/A1PK _ DEL-AZ  SPECTRUM  PARAMETERS,  LEFT  HALF 

C  IA2,  A2PK.  .  .  .  DEL-AZ  SPECTRUM  PARAMETERS,  RIGHT  HALF 

C  IE,  EPK. . DEL-EL  SPECTRUM  PARAMETERS 

C  PSL . SIDELOBE  LEVEL  (SEE  BELOW) 

C 


C  THE  REAL-TIME  OUTPUT  IS  THROUGH  /C5/ . 

C 

C  PEFF  ■  EFFECTIVE  POWER  RADIATED  FROM  RFSS  ARRAY  (FROM  SUB.  CLUTTR ) 

C  XR1  -  REAL  ARRAY  CONTAINING  SUM  CHANNEL  MODULATION  SIGNAL 

C  XII  *  IMAG  ARRAY  CONTAINING  SUM  CHANNEL  MODULATION  SIGNAL 

C  XR2  *  REAL  ARRAY  CONTAINING  DAZ  CHANNEL  MODULATION  SIGNAL 

C  XI2  *  IMAG  ARRAY  CONTAINING  DAZ  CHANNEL  MODULATION  SIGNAL 

C  XR3  *  REAL  ARRAY  CONTAINING  DEL  CHANNEL  MODULATION  SIGNAL 

C  XI3  =  IMAG  ARRAY  CONTAINING  DEL  CHANNEL  MODULATION  SIGNAL 

C 

C  IROT  IS  A  ROTATION  FLAG.  IF  IROT.  QT.  0  THE  AZ  AND  EL  ARRAYS  ARE 
C  INTERCHANGED. 

C 

C  THIS  PROGRAM  GENERATES  A  CONSTANT-LEVEL  SIDELOBE  SHOULDER  OF  (SPK-PSL) 
C  DB  BELOW  THE  SUM  SPECTRUM  PEAK. 

C 

C  ARRAYS  S,  A,  AND  E  ARE  FOR  TEMPORARY  STORAGE  OF  THE  SUM,  AZ-  AND  EL- 
C  DIFFERENCE  SPECTRA.  THEY  MUST  BE  DIMENSIONED  AS  LARGE  AS  2*NFFT. 

C 

C  THE  DIMENSIONED  ARRAYS  IN  COMMON  /C5/  MUST  ACCOMMODATE  NFFT+NOVLP 
C  WORDS  (SEE  TIMESQ). 

C 

COMMON  /C3/  NFFT 

COMMON  /C4/  IS(6)>  IA1 (6),  IA2(6)>  IE(6)«  SPK,  A1PK,  A2PK,  EPK,  PSL 
COMMON  /C5/  PEFF. 

I  XR1 (320),  XII (320) ,  XR2<320>, XI2(320),  XR3(320>,  XI3(320) 

COMMON  /C7/  S( 512) , A( 512) , E ( 512) 

COMMON  /LIM/  NFFT2 
,  NFFT2»NFFT*2 
C 

C  INITIALIZE  ARRAYS  TO  SIDELOBE  LEVEL 
C 

CALL  XMIT(-NFFT2, PSL, S> 

CALL  XMIT(-NFFT2» PSL,  A) 

CALL  XMIT(-NFFT2, PSL, E) 
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GENERATE  LOG-AMPLITUDE  SPECTRA 

CALL  CURVE<S, IS, SPK) 

CALL  CURVE<A, IA1, A1PK) 

CALL  CURVE< A,  I  A2> A2PK > 

CALL  CURVE < E,  IE,  EPK) 

CONVERT  TO  LINEAR  AMPLITUDE  AND  FOLDOVER 
SN0RM=*0. 

CALL  CONVRT  <  S,  SNORM ) 

CALL  CONVRT (A,  SNORM) 

CALL  CONVRT<E,  SNORM) 

GENERATE  TIME  SEQUENCES 

CALL  TIMESG< S, XRl, XII) 

IF <  IROT.  GT.  0)  CO  TO  35 
CALL  TIMESQ(A, XR2, XI2) 

CALL  TIMESQ(E. XR3,  XI3) 

CO  TO  40 

35  CALL  TIMESQtA,  XR3,  XI3) 

CALL  TIMESQ<E,  XR2,  XI2) 

40  RETURN 
END 
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SUBROUTINE  CURVE(S, II.  PK> 


C  THIS  SUBROUTINE  GENERATES  A  CURVE  IN  ARRAY  S  BETWEEN  SAMPLES  I I < 1 > 
C  11(5).  THE  CURVE  IS  DEFINED  IN  FOUR  SEGMENTS  BETWEEN  FIVE  POINTS 

C  WHERE  THE  POINTS  ARE  GIVEN  BY . 

C 


c 

SAMPLE 

VALUE 

c 

c 

II  <  1 ) 

0. 

c 

11(2) 

PK/2. 

c 

11(3) 

PK 

c 

11(4) 

PK/2. 

c 

11(5) 

0. 

c 

C  II(6)-ITYP 

IS  A  PARAMETER 

THAT  DEFINES  THE  TYPE  OF  SEGMENTS . 

C 

C 

ITYP 

SEGMENTS 

C 

C 

0 

LINE> LINE,  LINE/  LINE 

C 

1 

LINE/ PARAB,  LINE,  LINE 

C 

2 

LINE/ PARAB,  PARAB,  LINE 

C  THE  PARABOLA  IS  ALWAYS  CONSTRAINED  TO  HAVE  ZERO  SLOPE  AT  THE  CENTER 
C  POINT. 

C 

DIMENSION  S(l),  11(6) 

COMMON  /ABC/  A.  B,  C 


C 

c 

c 


c 

c 

c 


IF(PK.  LE.  0.  )  RETURN 
P2-  5*PK 
ITYP«II(6) 

COMPUTE  FIRST  SEGMENT 

c-o. 

S«P2/AMAX0( 11(2)— II(1)> 1) 

A— B*I  1(1) 

CALL  PARAB (5* 1 1 ( 1 ) *  11(2)) 

COMPUTE  SECOND  SEGMENT 

IF ( ITYP.  GT.  0)  GO  TO  20 
C-0. 

8»P2/AMAX0(II(3)-II(2).  1) 
A«P2-B*II<2) 

00  TO  25 

20  C— P2/AMAX0<  II  <3>-II(2>,  1)**2 
B— 2.  *C*II(3> 

A»P2-B#T I (2) -C*I I (2)**2 
25  CALL  PARAB (S, II(2)+1* 11(3)) 


AND 


O  O  O  CiCiCi 


COMPUTE  THIRD  SEGMENT 

IF ( ITYP.  GT.  1 )  GO  TO  30 
C*0. 

B*-P2/AMAX0( 11(4)— II(3)j  1) 
A*PK-B*II (3) 

GO  TO  35 

30  C»-P2/AMAX0(I1(4)-II(3),  1)**2 
B— 2.  *C*I  I  <  3 ) 
A»PK-B*II(3)-C*II(3>**2 
35  CALL  P ARAB (S»  II (3) +1.  II (4)  ) 

COMPUTE  FOURTH  SEGMENT 

C«0. 

B— P2/AMAX0(II(5)-II<4>,  1) 
A»P2-B*II (4) 

CALL  PARAB(S,  II(4>+1,  11(5) ) 


RETURN 
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SUBROUTINE  PARAB  <S<  III*  1 12) 

THIS  SUBROUTINE  GENERATES  A  PARABOLA  OF  THE  FORM. 

SID  *  A  +  B*I  +  C*I**2 

DIMENSION  S<1) 

COMMON  /ABC/  A,  B,  C 
COMMON  /LIM/  LIM 
I l=MAXO ( III.  1) 

I2=MINO( I 12* LIM) 

IFIIl.GT.  12)  RETURN 
DO  20  1  =  11,  12 
S(I)-(C*I+B)*I+A 
20  CONTINUE 
RETURN 
END 


SUBROUTINE  CONVRT ( S#  SNQRM ) 

C 

C  IN  THIS  SUBROUTINE  WE  CONVERT  THE  SPECTRUM  IN  ARRAY  S  OF  LOG-AMPLITUDE 
C  UNITS  TO  LINEAR-AMPLITUDE  UNITS.  IN  ADDITION,  TWO  PRF  INTERVALS  ARE 
C  FOLDED  (SUMMED).  THE  INPUT  S-ARRAY  IS  OF  LENGTH  2*NFFT  AND  THE  OUTPUT 
C  S-ARRAY  IS  OF  LENGTH  NFFT.  AN  APPROXIMATION  IS  USED  IN  THE  D0-20  LOOP 
C  FOR  A  FASTER  COMPUTATION. 

C 

C  FOR  THE  SUM  SPECTRUM  SET  SNORM=ZERO  ON  INPUT  (DO  NOT  USE  THE  CONSTANT 
C  0.  ,  HOWEVER).  THE  NORMALIZATION  FACTOR  SNORM  WILL  BE  RETURNED.  FOR 
C  THE  DIFFERENCE  SPECTRA  USE  THIS  SAME  VALUE  ON  INPUT. 

C 

DIMENSION  S<1) 

COMMON  /C3/  NFFT 
DO  20  K=1,NFFT 
S(K)=EXP(S(K)+S(K+NFFT) ) 

20  CONTINUE 

IF  (SNORM.  GT.  0.  )  GO  TO  35 
SUM-O. 

DO  30  K=l»  NFFT 
SUM=SUM+S ( K ) **2 
30  CONTINUE 

SN0RM=1.  /SORT (SUM) 

35  DO  40  K=*l»  NFFT 
S(K)«S(K)*SNORM 
40  CONTINUE 
RETURN 
END 
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SUBROUTINE  TIMESQ<SA.  XR,  XI) 

IN  THIS  SUBROUTINE  WE  GENERATE  A  CORRELATED  TIME  SEQUENCE  THAT  HAS  AN 
AMPLITUDE  SPECTRAL  FUNCTION  DESCRIBED  BY  ARRAY  SA  OF  LENGTH  NFFT.  THE 
METHOD  USED  IS  THE  ON-LINE  FFT  APPROACH  DESCRIBED  IN  SECTION  6  OF  MR  I 
REPORT  132-44,  WHERE  SUCCESSIVE  CALLS  TO  TIMESQ  WILL  CREATE  CONTIGUOUS 
TIME  SEQUENCES  BASED  ON  OVERLAPPING  FFTS  IN  THE  FREQUENCY  DOMAIN. 

ON  INPUT . 

NFFT  =  SIZE  OF  FFT  <=2** INTEGER) 

NOVLP  -  NUMBER  OF  SAMPLES  THAT  THE  FFTS  OVERLAP  (NFFT/B  IS  OK) 

SA(K)  =  SPECTRAL  AMPLITUDE  OF  KTH  SAMPLE,  K=1,...,NFFT 

ON  OUTPUT . 

XR(K)  =  KTH  SAMPLE  OF  IN-PHASE  SIGNAL,  K=1,...,NF 
XI (K)  -  KTH  SAMPLE  OF  QUADRATURE  SIGNAL,  K=1,...,NF 

WHERE  NF=NFFT— NOVLP. 

THE  ARRAYS  MUST  3E  DIMENSIONED  AS  LARGE  AS . 


SA 

NFFT 

XR 

NFFT +NOVLP- 1 

XI 

NFFT +NOVLP— 1 

A 

NOVLP- 1 

NOTE  THAT  THE  INFORMATION  IS  STORED  IN  THE  XR  AND  XI  ARRAYS  FROM 
NFFT+1  TO  NFFT+NOVLP+1  THAT  WILL  BE  USED  ON  THE  NEXT  CALL.  PRIOR  TO 
THE  FIRST  CALL  TO  TIMESQ  THE  INITIALIZATION  SUBROUTINE  (TIMSQI)  MUST 
BE  CALLED. 

THE  WORK  ARRAYS  WRKR  AND  WRKI  NEED  BE  DIMENSIONED  ONLY  IF  NFFT.  NE. 
2<-*  INTEGER,  IN  WHICH  CASE  THE  DIMENSION' SIZE  MUST  BE  NFFT. 

DIMENSION  WRKR < 160),  WRKI < 160) 

DIMENSION  SA( 1 ) , XR ( 1 ) ,  XI ( 1 ) 

COMMON  /C3/  NFFT, NOVLP,  NF,  A<1) 

CALL  GAUSS4 ( SA, XR, XI,  NFFT) 

•  CALL  FOURTtXR,  XI,  NFFT,  1,  0,  1,  WRKR,  WRKI) 

N1 -NOVLP- 1 

DO  20  K-1.N1 

,  XR(K)-A<K)*<  XR(K)-XR< K+NFFT) >+XR< K+NFFT) 

XI(K)-A(K)*(XI<K)-XI (K+NFFT) )+XI (K+NFFT) 

XR < K+NFFT )«XR(K+NFFT-NOVLP) 

•  X I ( K+NFFT ) -X I <  K+NFFT-NOVLP ) 

20  CONTINUE 

RETURN 

END 


A 
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+  SUBROUTINE  GAUSS4< A,  XR, XI,  N> 

C 

c  GENERATES  RANDOM  PHASOR  COMPONENTS  (XR(K)>  XI <K> )  OF  ZERO  MEAN  WITH  THE 
C  AVERAGE  POWER  OF  XR<K)**2+XI <K)**2  GIVEN  BY  A(K)**2  FOR 
C 

'  c  THE  RANDOM  NUMBERS  ARE  ACCESSED  FROM  RECIRCULATING  TABLES.  NTR  AND 

:•  C  NTI  ARE  THE  CIRCULATION  PARAMETERS  AND  NMAX=MAX(N).  WE  SHOULD  CHOOSE 

C  NTR  AND  NTI  RELATIVELY  PRIME. 

C 

C  ARRAYS  TR  AND  TI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTR+NMAX  AND  NTI+NMAX 
C  RESPECTIVELY. 

(ft  C 

C  IF  IRSET* 1  IN  COMMON  / GS /,  THE  POINTERS  WILL  BE  RESTORED  ON  EXIT  TO 
C  THEIR  ORIGINAL  VALUES  WHEN  GAUSS4  WAS  ENTERED.  THIS  MEANS  THAT  THE 
C  SAME  GAUSSIAN  SEQUENCE  (EXCEPT  FOR  THE  WEIGHTS  A)  WILL  BE  GENERATED  ON 
C  THE  NEXT  CALL. 

C 

•  C  PRIOR  TO  CALLING  GAUSS4#  THE  RANDOM  NUMBERS  IN  ARRAYS  TR  AND  TI  MUST 
C  BE  INITIALIZED  WITH  A  CALL  TO  GAUSSI. 

DIMENSION  A<1), XR<1), XI ( 1 > 

COMMON  /GS/  IRSET,  NTR,  NTI,  NMAX. KR, KI ,  TR < 1253),  TI< 1353) 

IF(N.  GT.  NMAX  )  STOP 
DO  35  1*1.  N 
XRU)*A<I)*TR(XR) 

XI ( I )=A( I )*TI (KI ) 

KR-KR+1 
KI-KI+l 
35  CONTINUE 

IF(  IRSET.  NE.  1)  GO  TO  40 
KR*KR-N 
KI-KI-N 

40  IF<KR.  GT.  NTR)  KR=KR-NTR 
IFCKI.  GT.  NTI)  KI-KI-NTI 
RETURN 
END 


! 


